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The calculation of the elastic scattering properties of mixtures composed by rubidium and potas- 
sium atoms is reported and compared with experimental results in detail. The improved potentials 
for both molecular states the singlet x 1 ^ + and the triplet a 3 J^ + of the RbK are presented, and the 
scattering lengths a t and the effective range r e are calculated with the help of this potential using 
WKB and Numerov methods for Rb-K in the triplet and singlet state. In addition, the convergence 
of these scattering properties as the depending on a Kq parameter is investigated using Quantum 
Defect Theory. The evaporative cooling other results that include the cross section and the spin- 
charge cross section, the rate coefficient as a function of the energy and the ultra-low temperature 
are also presented in this study. 
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I. INTRODUCTION 

Collisions of atoms at ultracold temperatures have re- 
ceived considerable attention because of their importance 
in the cooling and trapping of atoms Qand molecules 
and Bose Einstein condensation by Sympathetic cool- 
ing Ultracold mixtures of alkali metal atoms have 
recently been the subject of intensive experimental re- 
search, therefore the interest in the detailed knowledge 
of ultracold mixtures of alkali atoms has theoretically in- 
creased 0,H[- Even, while collisions either between alkali 
often same species or between different isotopes have of- 
ten been well characterized , collision between different 
alkali species have a still intensive exploring field. The 
condensed atoms and molecules may form a two-species 
Bose-Einstein condensate in which a new class of coher- 
ent phenomena is predicted to occur [2j. 

In an ultracold mixture, the interaction between the 
components is described in terms of a very important 
parameter, the interspecies s-wave scattering length. 
The character of the interaction between different alkali 
species in low-temperatures is determined by signs (±) 
and magnitude of the scattering length a: for bosonic 
atom with a>0 Bose condensate is stable, for a<0 it is 
unstable p|. This parameter also determines the effi- 
ciency of the sympathetic cooling procedure Q ■ In addi- 
tion, knowledge of the scattering length is important in 
order to predict macroscopic properties of these mixed 
systems, such as the stability and phase separation of bi- 
nary degenerate systems, or microscopic properties, such 
as the occurrence of resonances and the cross section and 
the spin-charge cross section in the interspecies collisions. 

Accurate calculations of s-wave scattering length and 
effective range for diatomic potentials are important due 



'Electronic address: 
t Electronic address: 



ozturkm® 
suleymant! 



jazi.edu.tr 
igazi.edu 



to the elastic collisions in these temperatures dominated 
s-wave scattering. In recent years, these quantities have 
been calculated by several theoretical methods 0, 0, 0, 
Il2l Il3l ] . The accurate determination of scattering length 
depends on choosing the interaction potentials. Collision 
processes at near- zero temperatures (~ fiK) are sensitive 
to the details of the interaction potentials between the 
colliding systems. 

Previous theoretical and experimental calculations of 
the scattering length for Rb+K interaction have ap- 
peared in the literatures 0, The purpose of this 
paper, we are to calculate scattering lengths, effective 
ranges and phase shifts-cross section for the collisions of 
Rb+K at ultra-cold temperatures with the Semiclassical 
method (WKB) 01 an d also solve the Schrodinger equa- 
tion using Numerov method numerically |l6j |. We have 
investigated the convergence of some scattering proper- 
ties obtained with help of these methods as the depen- 
dence of a K° parameter using Quantum Defect The- 
ory The interaction of rubidium-potassium atoms 

in states ( l ^2g) ano ^ ( 3 J2u) * s described by short and 
long range potentials which includes the cutoff functions 
introduced to truncate the unphysical short-range contri- 
bution of the polarization potential near the origin[r?j|. 
Cutoff functions depending on the cutoff radius are ad- 
justed by comparing the expressed potential form to the 
curve calculated with ab-initio by Rousseau et.al. |18| . 
Some results we have found are compared with theoret- 
ical values in the ref. 0] and the experimental result 
0. 

This paper is designed as the follow form, we discuss 
the computational process in Sec. 2. The interaction 
potentials, which play a crucial role in evaporate cooling 
and also determine many properties of the condensate, 
are given in Sec. 3. The ultracold scattering proper- 
ties(the scattering length, effective range, phase shifts 
and the rate coefficients) are computed for collisions of 
the Rb-K isotop systems in the singlet and triplet states 
at ultra-cold temperatures in Sec. 4 and 5, respectively, 
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and the convergences of the scattering range and the ef- 
fective range are also investigated in detail. We conclude 
with brief conclusions. 

II. PROCESSES 

The scattering length is defined from the asymptotic 
behavior of the solution of the radial Schrodinger equa- 
tion at zero energy: 



value 6 for interaction of mixing or same species atoms 
[T^ | . A similar calculation of the semiclassical phase at 
negative energies, E = — l/(2^a 2 ) yields <p(E) = <p — 
v /7F7 2 /"a-( n - 2 )/ T T(i+l/n)/[(n-2)r(l+l/n)].It is often 
used for studying quantization of vibrational levels near 
the dissociation limit in diatomic molecules pcj. 

The effective range r e can be written in terms of zero- 
energy solutions of the partial- wave equation if yo{r) is 
the solution of the partial-wave equation at k = for 
zero-potential, and normalized as 



2Sab + 1 



1(1 + 1) 



Vi(r) = o, 



(1) 

in atomic unit. Here ji is the reduced mass of the atomic 
system and V r (r) 2lSAf3+1 is the molecular potential be- 
tween two atoms. Asymptotic behavior of the wave func- 
tion is 



Vo 
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as k — ► 0, and if yo( r ) is the normalized solution of the 
radial equation Q for the s-partial wave at zero energy 
(for non-zero potential) as y(r) ~ Voi r) at r — > oo. Then, 
effective range can be written as [la, |21| 



y (r) =s r + s 



(2) 



as r — > oo, where s and so are constants. The scattering 
length is given by 



so 



(3) 



these constants are obtained by G. F. Gribakin using et 
al. using WKB method from semiclassical behavior and 
with help of the exact solution at zero-energy [l^ . Eq. 
© is transformed to the form 



a {1 — tan{<^ — 7r/8}} 



(4) 



where a is the "mean" or "typical" scattering length de- 
termined by the asymptotic behavior of the potential 
through the parameter 7 = \/2fm (/U is reduced mass, 
and k — cq is Van der Waals constant) for atom-atom 
interactions: 



- r ( 3 / 4 ) PT 

a = f(i7ir 27 



(5) 



Scattering length also depends a semiclassical phase ip 
calculated at zero energy from classical turning point ro 
where yijo) = 0, to infinity, 
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{Vo( r ) - V {r)}dr. 



(8) 



This integral converges provided y(r) approaches j/o( r ) 
rapidly enough as r — > 00. This requires V(r) to de- 
crease faster than r -5 . For — cg/r™potential type, this ex- 
pression is adjusted by Flambaum et al.j^ using WKB 
method and the exact solution of the radial equation at 
zero-energy as 
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The mean effective range may be easily obtained by 
replacing eq.Q into eq. © as 
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The low-energy scattering is dominated by the contri- 
bution 1 = 0, At values of k close to zero, the I — phase 
shifts £n can be represented by a power series expansion 
in k 



kcotSn = -- + -r e k 2 + o (k 3 ) . 

a 2 v ; 

This expression is known as Bethe formula. 



III. INTERACTION POTENTIALS 



(11) 



The total number of vibrational levels with zero orbital 
angular momentum was derived from ip, Nm — {p/ft — 
(71 — l)/2(n — 2)} + 1, where n is are given by integer 



The interaction between the atoms has two features: 
First, the potential at large distances behave as an in- 
verse power of the interatomic distance, V(r) = —n/r n 
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with n = 6 for spherically symmetric two atoms in their 
ground state. The asymptotic parameter k = Cq is known 
quite well for most atomic pairs of interest. Second, 
for atoms other than hydrogen and helium the potential 
curve is usually quite deep and also the electron-exchange 
part of the atomic interaction is repulsive at smaller dis- 
tance as for the singlet and triplet terms of alkali- metal 
atoms. "Deep" here means that the wave functions of 
the atomic pair oscillate many times within the poten- 
tial well, even at very low collision energies. The inter- 
atomic potential supports a large number of vibration 
levels. This latter feature enables one to use the semi- 
classical approximation to describe the motion of atoms 
within potential well 

The mixing alkali-metal atoms have unpaired electron 
spin s = 1/2 in the ground state Therefore, the interac- 
tion between two different alkali atoms in the ground 
state results in the formation of a diatomic molecule, 
which is described by two terms corresponding to the 
total spin s a b = s a + s& of the system equals to 1 or 0. 
The state with s a b = and 1 corresponds to the sin- 
glet state and triplet state, respectively. The probability 
of the spin exchange is determined by both the depen- 
dence of potential curves on the internuclear distance r 
and the splitting between singlet and triplet terms. The 
exchange interaction is analytically calculated using the 
surface integral method by Smirnov and Chibisov |24j . 
which yields 

v exc .(r) = l(v u -v g )=r a e(- f3r h. (12) 

Where a, (3 and C are constants. In this equation, J is 
given by ^2 n f n Jn^ n , here, f n are fit constants adjusted 
by fit program for the best R square in this paper, and J n 
are the integral constants. For Rb-K, these values of con- 
stants taken from the paper of Smirnov and Chibisov are 
given in Tabled The long-range part of the interaction 
potential is given by 

Eeven C m 

this term is the sum of Van der Waals terms, multiplied 
by the relevant cutoff function f c (r) to cancel 1/r" di- 
vergence at small distances. The area under r axis of the 
interaction potential is very important to specify the ul- 
tracold scattering properties, and its value is calculated 
from the eq. JHJ). For any different two potential, when it 
has same semiclassical phase shifts, the scattering prop- 
erties aren't affected from the change of other parameters 
of potential such as the depth of potential [25j . It is the 
fact that theoretical potential should actually adopt with 
ab- initio one. Therefore, in the singlet state, the some 
sensitive Wan der Waals parameters taken from ref. [2|| 
was multiplied with fit parameter Cj . As shown in the 
Figffl it concludes that the potentials are excellent adapt 
due to R square, and we may also name the Cj as the cor- 
rected constant for Wan der Waals parameters in singlet 



and triplet state, and this technique may be useful for 
other studies. Cg , C& and Cio parameters were taken 
from A. Derevianko et al.'s paper (Ref. 26] ). The cut- 
off function is analogous to that used for the H-H 3 E+ 
potential [27]]. 

fc(r) = £ (r - r c ) + £(r c - r)e~^-^\ (14) 

Where £(z) is unit step function £(z) = 1(0), when r > 
(<)0. There, the only free parameter for the potential 
sum is the cutoff radius r c , governing the decrease of the 
function f c (r) for r < r c . The potential of the collisions 
Rb+K consists of above the expressed two terms , the 
sum of Eqs. 1|12|) and (|13fl . is given as 

V 2sAB+1 {r) = V(r) long - range + (~l) s ^V(r) exc . (15) 

The magnitude of the cutoff parameters r c has been 
adjusted by comparing this potential with the experi- 
mental curve in triplet state obtained by Rousseau et.al. 
[l8| using ab-initio for the large range. 

The most relevant potential curves for different r c val- 
ues are represented in Fig^ together with the experi- 
mental calculation and their R square values (inset). Ow- 
ing to R square value, there is in excellent agreement be- 
tween experimental and theoretical curves for the cutoff 
radius determined from the fit curve. 



IV. THE CONVERGENCE AND 
CALCULATION SCATTERING LENGTH AND 
EFFECTIVE RANGE OF Rb + K 

In this section, the scattering length and the effective 
range for the collisions of Rb+K in singlet and triplet 
state at ultra-cold temperature have been calculated with 
semiclassical method (WKB) ^] and also determined by 
solving the Schrodinger equation by Numerov method 
[l(| for the potential given with Ea. (|15|) . Also, we in- 
vestigated the convergence of these scattering properties 
as a function of energy parameter(i4T°) using Quantum 
Defect Theory [13. 

In third column of Table [D] the results of the scatter- 
ing length calculated from Eq. (£1J with help of the WKB 
method are represented for the collision of mixtures com- 
posed by rubidium and potassium atoms in singlet and 
triplet state described with the potentials in eq. (jT5|) char- 
acterized by the term of cutoff radius. Some parameters 
are important in the determination of a(a.u.) for WKB 
method. We also obtain these parameters, the mean scat- 
tering length a and the zero-energy semiclassical phases 
ip, given in Table ITTll using values of asymptotic parame- 
ter 7. Here, the mean scattering length with asymptotic 
behavior of the potential in ea. i|15|) has been calculated 
from equation JSJ. 

For calculation numerically scattering length , we ob- 
tained s-wave phase shift So using the potentials in 
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ea. (|15fl . then we do it by solving radial equation using 
Numerov algorithm at e = fc 2 /2/i in atomic unit and find- 
ing Sq from the asymptotic behavior of the wave function 
y(r) — sin(fcr + 5q). The phases at small momenta k are 
used to extract the scattering length numerically from 

a = — linifc ,o (d5o/dk) The scattering lengths obtained 

via this method are given in the fourth columns of Table 
ITT1 The results obtained in two methods are compared 
with values found in the literature . 

From the measured thermalization rate, using the 
model presented in Ref. 28], the large value for the zero- 
energy s-wave scattering length are obtained using con- 
ventional double magneto-optical trap (MOT) apparatus 
[2^, l30| . G. Modugno, G Ferrari and et. al have es- 
timated a direct measurement of the scattering lengths 
for the 41 K- 87 Rb pair 0. They also deduced the colli- 
sional behavior of other isotope combinations by the help 
of mass scaling in the experimental columns of Table [H] 
[3. Then, for the 40 K- 87 Rb and 41 K- 87 Rb pairs, Simoni 
et al. investigated the scattering properties in the pres- 
ence of several magnetic field |3lj . 

As shown in the Table HT1 the scattering lengths calcu- 
lated for the collisions of rubidium- potassium pairs is in 
excellently agreement with the experimental scattering 
length when V(r) potential has the excellent adopt with 
ab initio curve due to the best R square value. 

Even our results shows that the stability of large 
condensates requires repulsive interactions (positive o), 
whereas for attractive interactions (negative a) it is un- 
stable, only a finite number of atoms can be found in 
condensate state in a trap. As shown Tablo[n] the scat- 
tering lengths found as negative value due to changes in 
the V(r), lead to a condensate triplet state where the 
number of atoms is limited to a small critical value de- 
termined by the magnitude of a. In contrast, we have 
observed the positive scattering lengths that produce sta- 
ble condensates. Also, from the scattering length given 
as a function of the semiclassical phase shift in the eq. 
(@J, we say that the minor changes in phase shifts have 
strongly affected the scattering amplitude. It concludes 
that the ultra-cold scattering depends on the potential 
changes given by phase. It must be remarked that our 
emphasis for 41 K- 87 Rb pair is weakly dependent on the 
number Nm of bound states supported by ab-initio po- 
tential. The numbers of bound states Nbd have been com- 
puted in the intervals 108-110 and 32-33 for both singlet 
and triplet states, respectively. The numbers of bound 
states of the singlet state is consistent with the result 
(Nb = 32) obtained ref.0] In addition, the mass scalling 
from one isotope to the others depends more strongly on 
n s . The calculations in this paper shows that due to be- 
ing boson- fermion, of great interest also is the 40 K- 87 Rb 
pair, for which we find a negative and probably large in- 
terspecies scattering length agreed with other result in 
Table llTI which implies an efficient sympathetic cooling of 
the fermionic species down to the degenerate regime [l4j . 

Effective ranges calculated using Numerov and WKB 
methods in the function of phase for the triplet rubidium- 



Potassium potentials are presented in Table llVl together 
the calculation of semiclassical phase shifts. The close 
agreement between the calculations of r e from Eqs.JHJ 
and (|llfl confirms the accuracy of the numerical integra- 
tion of the partial- wave equation. The size of the scatter- 
ing lengths and effective ranges is closely related to the 
position of the last vibration bound states of the energy 
curves, as can be anticipation by inspection of Eq.@ 
and number of vibrational levels n s (tp/ir) which, consis- 
tent with Levinson's theorem, show that as the binding 
energy of the highest level tends to zero, the scattering 
length tends to zbinfinity |l5|. 

Analytical calculations of scattering lengths are im- 
portant in investigation of convergence of the scattering 
length and the effective range. In many approaches used 
to solve collision problems in atomic physics, the tree di- 
mensional configuration space is divided into two regions 
separated by a spherical shell (a core boundary) of radius 
p. In the inner region (r< p) the short-range interaction 
between two colliding particles is very complicated and 
a scattering equation must be solved independently for 
each combination of particles. In contrast, it p is cho- 
sen sufficiently large the scattering problems in the outer 
rcgion(r> p) may be reduced to potential scattering with 
the long-range potential accurately approximated by sim- 
ple analytical expression. A numerical solution in this 
region is usually easily approachable [13] ■ One of these 
methods is Quantum Defect Theory described by Gao 

Quantum defect theory of atomic collisions is pre- 
sented by Bo Gao. Based on the exact solutions of the 
Schrodinger equation for an attractive 1/r 6 potential, the 
theory provides a systematic interpretation of molecular 
bound states and atom-atom scattering properties and 
establishes the relationship between them. Bo gao finds 
a definition for the scattering length and the effective 
range and s-wave at zero energy as 



— 2ti- K u -l a 
ai=o - pr(i/4)]» -pn-Pe 

_ _ [r(i/4)l 2 M'+i q 

r '=0 - 3^ \K^~^ ) 



(16) 



where /3q = (2fiCe) 1 ^ 4 and K°is the analytic function 
of energy. 

Figure [21 show the scattering length and the effective 
range curve plotted with the help of eqs. (|16fl and val- 
ues (dotted points) obtained by the mean WKB formulas, 
eqs.© and (fT77f> . 

Our graphical representations show that as K° goes 
infinity, the QDT results converges excellently to values 
obtained by WKB method in Fig. [5] . Even, it may 
be said that the different calculating of the scattering 
properties in these two method gives similar results for 
any ultracold atom-atom collision (e.g. ref. 33]). 
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V. PHASE SHIFTS AND CROSS SECTION 

Scattering length describes behavior of the atom scat- 
tering at low energy, and phase shifts are important pa- 
rameter for it. Furthermore, at a low energy, E, of rela- 
tive motion the phase-shift is N^tt + 5, where Nm is the 
number of bound states. The phase shift can be calcu- 
lated by fitting the solution y(r) to sin(fcr) and cos(fcr), 
the asymptotic solutions of Eq. QJ. Let Si and Cj de- 
note sin(fcri) and cos(fcri), respectively. Fitting at and 



we find, for a small shift [ljjllj|, 



5 w tan((5) 



CiUj — Cjyi 



(17) 



where r^j = r + (i ,j)h (h is the step-length ). we show 
the numerical phase shifts (module tt) in Fig. [3]it is also 
known in more tex book that the low energy scattering 
is always dominated by the I = partial wave the corre- 
sponding phase shift being given by 



-ka 



(18) 



which is shown for all Rb+K collisions in Fig. [3J Here, 
to plot phase shifts we used the semiclassical scattering 
lengths. As momenta goes zero, for I = 0, the s-phase 
shifts are very close in the singlet state, it concludes 
that when a gas mixture composed of fermion-boson or 
boson-boson is cooled below a critical temperature, the 
mixing atoms condensates in the lowest quantum state. 
Moreover, if the atoms are bosons, a cloud of atoms oc- 
cupies the same quantum state at low temperature. If 
the atoms are fermions, cooling gradually brings the gas 
closer to being a 'Fermi sea' in which exactly one atom 
occupies each low-energy state [34j . 

As shown in fig. |3 their scattering lengths are unstable 
against collapse except the 40 K- 85 Rb. For these mixtures 
of composed by fermion-boson, the condensates are more 
attractive than that of boson-boson as reported in Table 
IIIII In this case, the fermionic cloud cannot shrink below 
a certain size determined by the Pauli Exclusion Principle 
due to bosonic cloudy 

Elastic singlet and triplet cross sections between two 
atoms, can be defined by 



S,T 



4tt 



E 

(=0 



(21 + 1) sin 2 <J 



2 X S,T 



(19) 



where S and T stand for singlet and triplet, respectively. 
We have obtained the phase shifts (Si) and cross sections 
(dai) for the collision 85 Rb+ 40 K with I = 2,4 values of 
the angular momentum and different energies using Nu- 
merov method. They are listed in Table ITvl this table 
shows that the Si and do\ for very low energy have large 
values lead to the being of more collisions. In the ultra- 
cold gas of the mixing alkali metal atoms dominated by 



the I = contribution, the cross sections are given by 

An a| T for boson-boson sys- 



S,T 



8tt ag rp and a e[ 



S,T 



terns and boson- fermion systems, respectively 35]. The 
spin-change cross sections are calculated from the singlet 
and triplet phase shifts by 



i=o 



l)sin 2 (^ 



6?) 



(20) 



In the limit of low energies, Eq. 1201 is transformed to 
a sc = 7r (oj, — a,j>) 10]. In addition, this parameter is 
known as spin-flip collisions lead to trap losses [3j|. The 
zero-energy cross sections for the selected mixing atoms 
are 3.29010 +11 ag(=9.21310~ 10 m 2 ) and 3.06510 10 ag 
(=8.58410~ n m 2 ) for singlet and triplet states of the 
40 K- 85 Rb, and 4.700 10 +11 a 2 (=1.316 10-° 9 m 2 ) and 
8.152 10+ 10 (=2.28310- 10 m 2 ) for singlet and triplet 
states of the 39 K- 87 Rb, respectively. Spin-charge cross 
sections are also 7.45910 +10 a 2 l (=2.08710- 10 m 2 ) in the 
40 K- 85 Rb and 2.27610 +10 a 2 1 (=6.37310- n m 2 ) in the 
39 K- 87 Rb. 

We assume the velocity distribution is Maxwellian 
characterized by a kinetic temperature T , and we de- 
fine mean elastic and spin-change cross sections by 



oo 

d{T) = Jkkf I a w Eei ~ E/kBT)dE - ( 21 ) 



Assuming equal temperatures of the each pair of iso- 
topes in the Rb-K vapor, the corresponding rate coeffi- 
cients are given by 



S,T,SC 



TTfJ, 



a(T). 



(22) 



We show values of the cross sections as depending on 
the energy in the (m 2 ) in Fig. Q]for E up to 10 _2 (a.u.). 
and list values of the spin-change 2s ab + 1 in cm 3 s _1 in 
the Table IVll for T up to IK. 

Figure 4 gives cross sections for some interaction of 
atoms. They show that in the low energies, less partial- 
wave contributes to scattering . As shown in figure 
fermion-boson 85 ' 87 Rb+ 40 K collisions display different 
properties, it may be said that 85 Rb+ 40 K is more stabil- 
ity than 87 Rb+ 40 K, because later has less partial waves 
with decreasing values of the energy for both singlet and 
triplet state and its total cross sections are in very low 
values in the same higher energy region. 



VI. CONCLUSIONS 

The elastic scattering properties for the collision of 
the two mixing atoms in a limited range of moment 
(k<a _1 wO.Ola.u. for Rubidium-Potassium isotopes) at 
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low temperatures are sensitive to the details of interac- 
tion potentials. Scattering properties such as the scat- 
tering length, the effective range, the phase shift and the 
cross section have been computed using semiclassical and 
a numeric methods for the V(r) potentials as dependence 
of cutoff radius adjusted by comparing with the exper- 
imental potential for ultra-cold rubidium-potassium iso- 
topes collision. We investigated the convergence of these 
scattering properties as the depending of core radius and 



K° parameter using Quantum Defect Theory. Cross sec- 
tion was obtained as ~1.010 _10 m 2 at low energy. The 
phase shifts, an intermission parameter for scattering 
length and effective range, has linear manner at small 
momenta k < a -1 . So, it has been determined by com- 
puting numerically and analytically for V(r) potential. In 
bigger momenta than the inverse mean scattering length, 
the phase shifts display parallel behavior as the depen- 
dence of semiclassical phase shift, So(k) ~ ip — ka. 
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TABLE I: Potential coefficients for Rb-K atoms in singlet and 
triplet state. 







3 + 


l(T 3 /i 


345 


-345 


10" 3 /2 


1.0 


-1.0 


10~ 3 /3 


1.0 


-1.0 


10~ 3 Ji 


5.0 


5.0 


lO" 4 ^ 


4.0 


4.0 




7.0 


7.0 


a 1 


5.256 


5.256 




1.119 


1.119 


io- 4 cf 


1.502 


4.274 


10" 5 C^ 4 


5.922 


5.922 


10- 7 C? 


6.726 


6.726 



1 Ref.[24], 2 Ref.[26j. 

3 The fit constant defined in the text isC;f = 3.514 
4 The fit constant defined in the text is Cg' T = 0.865 



8 



TABLE II: Scattering lengths a (in atomic units) for ultracold 
rubidium-potassium mixtures in singlet and triplet states. 





39 K 






40 R 






41 R 






o" 


a" 1 


a 1 


o" 




a 1 


a" 


a" 1 



Rb S 136.73 136.76 _ 161.84 161.807 _ 212.56 212.55 _ 



T 56.983 56.952 58_ 6 -49.364 -49.390 -38L 7 i 7 271.568 306.016 329^ 

87 Rb S -3.121 -3.197 _ 23.927 23.926 _ 54.459 54.456 _ 

T 36.854 36.894 3lL 6 6 -200.05 -200.04 -26ll_ 7 1 ° 59 167.29 167.28 163^ 
'Calculated using semiclassical method for the V(r) potential. 
"Calculated using Numerov method for the V(r) potential. 

*"Ref.[14], scattering lengths have been measured with help of the conventional 
cross-dimensional thermalization technique. 
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TABLE III: Mean scattering lengths a (in atomic units) and 
scmiclassical phase shifts for ultracold rubidium-potassium 
mixtures in singlet and triplet states. 

39 K 40 K 41 K 

85 Rb a s ,ip 93.65 339.254 94.056 342.202 94.449 345.072 

a t ,ip 68.393 101.089 68.690 101.967 68.977 102.823 
87 Rb a s ,<p 93.82 340.486 94.230 343.467 94.627 346.369 

at,ip 68.517 101.357 68.817 102.244 69.107 103.108 



TABLE IV: Effective ranges r e (in atomic units) for rubidium- 
potassium atom scattering. 



39 


K 


40 


K 


41 


K 


r\ 


r" 




r" 


r\ 


r" 


85 Rb S 155.323 
T 295.498 


155.321 
295.400 


140.841 
1534.41 


140.841 
1534.38 


139.505 
124.996 


139.502 
124.989 



87 Rb S 511.371 511.387 6637.97 6635.634 983.868 983.527 

T 838.647 838.413 386.47 386.535 103.87 102.995 
'Calculated using semiclassical method for V(r). 
"Calculated using Numerov method for V(r). 



TABLE V: Cross sections (in the m 2 unit ) and phase 
shifts calculated for V(r) potential using numerical Numerov 
method. 

"Rb-^K 



Lgio(E) L Factor 



85-39 85-40 85-41 87-39 87-40 87-41 



8 -14.99 2 10-* 5.17 4.93 5.86 3.28 4.48 2.69 

4 10~ 3 5.50 5.57 6.61 3.64 5.36 3.36 

-11.57 2 - 0.256 0.245 0.289 0.165 0.223 0.136 

4 - 0.271 3.261 0.322 1.562 0.265 0.169 

da -14.99 2 10~ 12 5.00 4.92 7.67 2.24 4.73 1.58 

4 10~ 12 3.63 3.55 5.61 1.62 3.43 1.14 

-11.08 2 10~ 12 4.57 4.52 6.63 2.15 4.33 1.53 

4 10~ 12 3.32 3.26 4.85 1.56 3.14 1.12 



TABLE VI: Singlet, Triplet and Spin-change rate coefficients 
for different isotops of the mixed Rb-K. The numbers in brack- 
ets denote multiplicative powers, a) 85 Rb- 39 K, b) 85 Rb- 40 K, 
c) 85 Rb- 41 K, d) 87 Rb- 39 K, e) 87 Rb- 40 K, f) 87 Rb- 41 K. 
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FIG. 1: Potential energy (dotted line) curves obtained as 
compared with experimental (solid line) potential curve for 
the best cutoff radii in the singlet and triplet states of the 
moleculer KRb. The insets show R 2 values and cutoff radii 
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FIG. 2: Convergence of the scattering length and the effective 
range calculated with QDT as a function of K parameter 
when K goes infinity. 
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FIG. 3: s-wave phase shifts calculated for the V(r) potential 
in the singlet and triplet states (Solid lines). Dotted lines 
show the curves calculated with 8q— -a k. 




FIG. 4: Cross-sections as a function of the energy in the spin 
charge, the singlet and triplet states of the mixing of atoms. 



